##Exercicio 5###

##5.1:Crie um grafico de dispersao entre dap1) e ht2) 

setwd("C:/Users/GUILHERME/Desktop/Vivian/Linguagem R/Aula5/Exercicios")
esaligna<-read.csv2("esaligna.txt",header=TRUE,sep=",",as.is=TRUE)
esaligna
str(esaligna)
class(esaligna)
head(esaligna)

1. legenda nos eixos com os nomes das variveis e suas unidades
plot(esaligna$dap~esaligna$ht,xlab="DAP(cm)",ylab="Altura(m)")

2.Marcacoes dos eixos para dentro da rea do grfico
plot(esaligna$dap~esaligna$ht,xlab="DAP(cm)",ylab="Altura(m)",tcl=0.3)

3.Apenas 2 eixos (formato "L")
plot(esaligna$dap~esaligna$ht,xlab="DAP(cm)",ylab="Altura(m)",bty="l",tcl=0.3)

4 e 5. ttulo informativo e tamanho das fontes maiores que o padro
plot(esaligna$dap,esaligna$ht,xlab="DaP (cm)",ylab="Altura (m)", main="DAP x ALTURA",bty="L",tcl=0.3,cex.main=1.5,cex.axis=1.5) 


##5.2. Dois grficos juntos
1. Use as variveis dap e talhao para construir dois grficos

setwd("C:/Users/GUILHERME/Desktop/Vivian/Linguagem R/Aula5/Exercicios")
esaligna<-read.table("esaligna.txt",header=TRUE, sep=",",as.is=TRUE)
esaligna
str(esaligna)
class(esaligna)

# Os grficos do DAP e talhao juntos
#Primeiro grfico
par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$talhao,xlab="Talho", ylab="DAP")
mtext("a")

#Segundo grfico
media.dap<-tapply(esaligna$dap,esaligna$talhao,FUN=mean)
media.dap
desvio.padrao.dap<-tapply(esaligna$dap,esaligna$talhao,FUN=sd)
desvio.padrao.dap
esaligna2<-data.frame(media.dap,desvio.padrao.dap)
esaligna2
var1<-esaligna2$media.dap+esaligna2$desvio.padrao.dap
var1
var2<-esaligna2$media-esaligna2$desvio.padrao.dap
var2
plot(unique(esaligna$talhao),esaligna2$media,xlab="Talho",ylab="DAP",ylim=c(min(var1), max(var2)))
segments(unique(esaligna$talhao),esaligna2$media,y1=var1)
segments(unique(esaligna$talhao),esaligna2$media,y1=var2)
points(unique(esaligna$talhao),var1,pch="---")
points(unique(esaligna$talhao),var2,pch="---")
mtext("b")


##Exercicio 5.3
arquivo<-read.table("arquivo.txt",header=TRUE, sep=",",as.is=TRUE)
arquivo
str(arquivo)
class(arquivo)
head(arquivo)
summary(arquivo)
par(mfrow=c(1,2))

#primeiro grfico
?plot
plot(y1~x1,data=arquivo,xlab="Log(Patch size)(ha)",ylab="Euclidean distances",tck=0.02,cex.lab=1,cex.axis=0.7,family="serif",pch=15,bty="L",xlim=c(0.5,3),ylim=c(0.5,3.0)) #teste1 
plot(y1~x1,data=arquivo,xlab="Log(Patch size)(ha)",ylab="Euclidean distances",tck=0.02,cex.lab=1,cex.axis=0.7,family="serif",pch=16,bty="L",xlim=c(0.5,4),ylim=c(0,3.0)) #teste2
plot(y1~x1,data=arquivo,xlab="Log(Patch size)(ha)",ylab="Euclidean distances",tck=0.01,cex.lab=1,cex.axis=0.6,family="serif",pch=17,bty="L",xlim=c(0.5,2),ylim=c(0,3.0)) #agora sim! 
linha<-lm(y1~x1,data=arquivo)
abline(linha)
mtext("a",family="serif",cex=1,at=2)


#Segundo grfico
par(bty="L",tcl=0.3)
boxplot(x2~y2,data=arquivo,ylim=c(0,3),names=c("Small","Medium Edge","Medium Interior","Large Edge","Large Interior","Control"),cex.main=1.2,cex.axis=0.5) #teste1:ficou com um outline
boxplot(x2~y2,data=arquivo,ylim=c(0,3),outline=FALSE,names=c("Small","Medium Edge","Medium Interior","Large Edge","Large Interior","Control"),cex.main=1,cex.axis=0.5) #agora sim!!!
mtext("b",at=6.5,cex=1,family="serif")
mtext("*",at=1,cex=1,line=-1,family="serif",font=1)
mtext("*",at=2,cex=1,line=-1,family="serif",font=1)
mtext("**",at=3,cex=1,line=-1,family="serif",font=1)
mtext("*",at=4,cex=1,line=-1,family="serif",font=1)
mtext("***",at=5,cex=1,line=-1,family="serif",font=1)